The Interplay of Spin and Charge Channels in Zero Dimensional Systems 
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We present a full fledged quantum mechanical treatment of the interplay between the charge 
and the spin zero-mode interactions in quantum dots. Quantum fluctuations of the spin-mode 
suppress the Coulomb blockade and give rise to non-monotonic behavior near this point. They 
also greatly enhance the dynamic spin susceptibility. Transverse fluctuations become important as 
one approaches the Stoner instability. The non-perturbative effects of zero-mode interaction are 
described in terms of charge ((7(f)) and spin (5(7(2)) gauge bosons. 
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The importance of electron-electron interactions is 
emphasized in low-dimensional conductors. In one- 
dimension interactions in the charge and spin chan- 
nels are separable Considering zero-dimensional quantum 
dots (QDs), the "Universal Hamiltonian" [3,0 scheme 
provides a framework to study the leading interaction 
modes: zero-mode interactions in the charge, spin (ex- 
change) and Cooper channels. While this Hamiltonian 
is simple, the physics involved is not at all trivial. The 
charge channel interaction leads to the phenomenon of 
the Coulomb blockade (CB). The exchange interaction 
leads to Stoner instability [3j, which, in mesoscopic sys- 
tems as opposed to bulk, is modified Attention has 
been given to the intriguing interplay between the charge 
and the spin channels. This is manifest, for example, in 
the suppression of certain Coulomb peaks due to "spin- 
blockade" 0. In a recent theoretical study the effect 
of the spin channel on Coulomb peaks has been analyzed 
employing a master equation in the classical limit. 

In this Letter we study both transport through a 
metallic grain and the dynamic magnetic susceptibility 
of the latter. Specifically we find that (i) the spin modes 
renormalize the CB, thus modifying the tunneling density 
of states (TDoS) of (hence the differential conductance 
through) the dot (cf. Fig.2 and Eq. I55)l . For an Ising- 
like spin anisotropy (represented by 1 — e) the longitudi- 
nal mode partially suppresses the CB. Transverse modes 
act qualitatively in the same way, but as one approaches 
the Stoner instability point (from the disordered phase), 
the effect of transverse fluctuations reverses its sign and 
acts towards suppressing the conductance (i.e., enhanc- 
ing the CB). This results in a non-monotonic behavior of 
the TDoS. (ii) The longitudinal magnetic susceptibility 
\ zz l|25(l diverges at the thermodynamic Stoner Instabil- 
ity point, while (I25fl is enhanced but remains finite. 
However, one notes that the static transverse susceptibil- 
ity is enhanced by the gauge fluctuations. 

Our study is the first full fledged quantum mechanical 
analysis of spin fluctuations and the charge-spin inter- 
play in zero dimensions. The non-perturbative effects 



of zero-mode charge interaction (e.g. zero-bias anomaly 
0) are described in terms of the propagation of gauge 
bosons ((7(1) gauge field) 0. Here we adopt similar- 
ideas to account for spin fluctuations described by the 
non-abelian SU(2) group. The Coulomb and longitudi- 
nal spin components are accounted for " exactly" , while 
transverse spin fluctuations are analyzed perturbatively 
(with the easy-axis anisotropy, e 0). These fluctuations 
become important as one approaches the Stoner insta- 
bility. Here we restrict ourselves to the Coulomb valley 
regime and ferromagnetic exchange interaction. 

Before proceeding we recall that beyond the thermo- 
dynamic Stoner instability point, 3th — A (A being the 
mean level spacing) the spontaneous magnetization is an 
extensive quantity. At smaller values of the exchange 
coupling, J C <J < Jtin finite magnetization shows up |l(, 
which, for finite systems, does not scale linearly with 
the size of the latter 0. Its non-self-averaging nature 
gives rise to strong sample-specific mesoscopic fluctua- 
tions. The incipient instability for finite systems is given 
by J c = A/(l + e) for an even number of s pins in the dot 
and J c = A/(l + e/2) for an odd number 

Hamiltonian and correlators. Our QD is taken to be 
in the metallic regime, with its internal dimensionlcss 
conductance g 1. Discarding both Cooper and spin- 
orbit interaction channels, the description of our metallic 
QD allows for only two other channels, namely charge 
and spin. While the charge interaction is invariant un- 
der (7(1) transformation, the spin interaction possesses 
a non-abelian 5(7(2) symmetry associated with the non- 
commutativity of the quantum spin components. 

The Universal Hamiltonian now assumes the form 



H — 22 e aal l (T a a:Cr + He + Hs 



(1) 



Here a denotes a single-particle orbital state with spin 
projection a. For simplicity, below we confine ourselves 
to the GUE case. The Hamiltonian He — E c (n — N ) 2 
accounts for the Coulomb blockade, E c = e 2 /2C is a 
charging energy, n the number operator; Nq stands for a 
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positive background charge tuned to a Coulomb valley. 
The Hamiltonian 



H 8 = -J 



E*« 



E^ + E^ 



represents the spins S aa ' =\ J2 a a a o& 



(2) 

- 1 ^ - -» a j interac- 
tion within the dot. Hereafter we assume strong easy 
axis anisotropy 8], e = J±/ J\\ < 1. In this case the spin 
rotation symmetry is reduced to SO(2). We will treat the 
terms of transverse and longitudinal (Ising) fluctuations 
independently. 

The Euclidian action for the model QJ is given by 



S = 



C{r)dT 



ch 



(3) 



Here ^ stand for Grassmann variables describing elec- 
trons in the dot. The imaginary time single particle 
Green's function (GF) is written as 



£[^# Q ^)Vw(T/)e s ^] 



where partition function Z(p) is given by 



Z(jjl) = / D[W>\ 



(4) 



(5) 



Employing a Hubbard-Stratonovich transformation 
with the bosonic fields (j> (for charge) and <E> (for spin) 



cxp f- J drHc^j = J D[4>] exp \ - J dr 
cxp(-J^ drHsj = J D[$]cxp(- J dr 



E n a ~ N 
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(G) 



we obtain a Lagrangian which includes a term quadratic 
in f = ^ Q ^cMq^q. Here we have used a spinor 
notation ^> a — (^-faV'J.a) & n d the matrix M Q is given by 



M„ = 



$ 2 



(7) 



Our goal here is to obtain the GF. We first add source 
fields to the Lagrangian 

£ AT = C + A* + *A + f I, 

and define the generating function Z as follows 

Z = Zifi)- 1 J D[**]D[^]exp (J%tUt(t) \ , (8) 

where Z{p) is a partition function JSJ) of the dot. The 
fcrmionic (2x2) and bosonic (3x3) matrix GFs are: 



d 2 Z 



dA^dA 



d 2 Z 



(9) 



with A — > 0, T — > 0. Here <J Q is given by Q aa = 
-(^^(t/)*^)) while V iU = -(T T *'*(r i )4 , '(r / )). 

Gauge transformation. We now apply a (non-unitary) 
transformation to gauge out both the Coulomb and 



the longitudinal part of the spin interaction M a = 
WM a W- x . We have ^=W(t)^ and ^=#W _1 ( r ) with 



W(t) = e 



- J6(t) 



e V(r) 

e-"( T ) 



(10) 



Here 9 (r/) accounts for the U(l) fluctuations of the charge 
(longitudinal) fluctuations, 



(<f>(r')~^)dr', n= (^(r')-^)dr'. 



(11) 



ac- 



In defining the gauge fields <po and $g one needs to . 
count for possible winding numbers (fc,m=0±l, ...) 

(3(f> a = / 4>{r)dT + 27rfc, = / $ z (r)dr + 2mm 

Jo Jo 

(12) 

In Ea. (|llf> initial conditions W(0) = 1 and periodic 
boundary conditions W(0) = W((3) are used. As a 
result, the diagonal part of the gauged inverse elec- 
tron's GF (M Q ) does not depend on the finite fre- 
quency components of fields. The off-diagonal part 
can be taken into account by a perturbative expansion 
in e < 1. We represent M a = (G^y 1 + with 



^ U| (r) 



(d T — £ a + i(f>o)l + $>QO- z and the self-energy 
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E$ = Ve ( < J ,_ e 2 ''cr+ + $+e- 2 ''cr-). We next calculate 
the Green's function 

Gf (r, p) = 5 aa , < SL 01 (r, fi) exp (ar, r + i9 T )) $(j> . (13) 

Hereafter ((•■■)tr) denotes Gaussian averaging over 



fluctuations (transverse fluctuations) of the bosonic field 
(<&,</>) and fL=^,+a^Q+i(j)Q. Integrating over all Grass- 
mann variables and expanding M a with respect to the 
transverse fluctuations, one obtains [l2j | 



c?$o^0o exp 



(%? (0o) 2 



TJ 4TE C 



JJd^# n exp^T^ 

n^O I n^O 



J 



4B„ 



^(Ti-r/.M + i^o+^g)-^ / dn... / dT 2k {(Q^) 2k ) tr QW{T 2k -T h ii + i<p + a^) 

k—1 



W- 1 {tj). 



(14) 




/3 



(17) 



The exchange interaction effectively modifies the charg- 
ing energy. For long-range interaction this correction is 
small E c / J^(fci?i) d_1 2] (X is a linear size of the d- 
dimensional confined electron gas, Uf is a Fermi momen- 
tum), while for contact interaction E c — J/4 = Q. 

Transverse fluctuations. The first non-vanishing dia- 
gram of our expansion l|14(l is depicted in Fig. lb. Here 



FIG. 1: First and second order Feynman diagrams con- 
tributing to electron's GF. Solid line represents Qaja-; double 
dashed line stands for combination of Coulomb and longitu- 
dinal bosons; single dashed line denotes a longitudinal boson 
while the zig-zag line represents ($ + (ti)$~(t2)}. 



where S1o(m) = —TlnZo, Zq is the partition function of 
the non-interacting electron gas. Also, computing the 
bosonic correlator (Eq J^J), we find 



1 



Z(jjl) 



x exp Tr log 



1 + Sl 01 ^ 



(15) 



In the spirit of 7] , the interaction of electrons with the 
finite-frequency charge and longitudinal modes (4> n , $*) 
may be interpreted in terms of a gauge boson 0] dressing 
the electron propagator (cf Fig. la). The exact electronic 
GF (depending on winding numbers i s given by 

G a An - r f ) = Q^An - r f , ^ s ^f\ (16) 

where the Coulomb- longitudinal U(l) gauge factor is 



5 ll( r )= 4T ^ 

n=iO 



E c - J/4 



sin 



^7 



(rf) 



S'SM = («{„ (1 - 6 T ) - (1 - n ( „, )0 T ) (18) 

the transverse correlator is considered in the Gaussian 
approximation 



(n)*-(7a)) = ^(n 



eJ 2 



T ^ + W^7Ty (19) 

In Ea . (|19fl the first term is a manifestation of the 
white noise fluctuations of the fields <& arising from the 
Gaussian weight factor (cf. Eq. l|15[lV The second term 
is related to the non-Gaussian factor in T>^ v and reflects 
the feed-back of (the ^-dependent) on T>. Note that 
the transverse components $ ± are always accompanied 
by the gauge factors e T2v , hence the longitudinal bosons 
contribute to the dynamics involving the transverse 
fluctuations. 

To proceed we now sum Ea. (114(1 over a. Perturbative 
corrections to the electron GF coming from transverse 
fluctuations are now expanded in e and summed up in 
the factor F± 

G(T)=Go(r)e' s ^F ± (T,e), (20) 



where 



F ± (r, £ ) = l + ^{/("V m) }- (21) 
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The effects of disorder are incorporated in the bare den- 
sity of states vo = 1/A. We denote r=r/— r*. At fi- 
nite temperatures we employ the conformal transforma- 
tion l/(Tj — Tf) — > 7r/(/3sin (7r(Tj — Tf)/ (3)). The factors 
{y( n )^( m )} refer to diagrams containing n + to 2? (zig- 
zag) line correlators, in n (to) of which we employ the 
first, <5-function (second, constant) term of Eq. I|19|l . The 

Here we calculate F± to 



,n+2m 



factor is 
order e 2 . 

In general, we may write the effective transverse gauge 
boson as D±(e,r) = exp (— S±(e,r)) Q. S± preserves 
the symmetry (in t) with respect to (3/2 to all orders of 
e. It can therefore be written as <Sj_(e, t — ^-). 

The first term in (|21|) (of order e) is given by 



(i) - 



£ 2 



(22) 



This contribution is of the same origin as that of the 
longitudinal boson part (namely it comes from the 
j <£> 2 d,T term in Eq. l[T§|l.l Along with the other 
to = terms in H21|l it can be exponentiated |14| . result- 
ing in J/4 — > J(l + 2e)/4 in the expression for S\\. For 
the isotropic model one obtains (1 -I- 2e)/4 — > S(S + 1). 

There are contributions to S± arising from the sec- 
ond term of Eci. (|19fl . As a result, the lowest, ~ e 2 
contribution (the {f^°'g^} term of the expansion l|21|l) 
leads to a non-Gaussian contribution to S± which is 
~ ~ 2[3(Js/ 1 e J) t 7 " — "^"] 2 ' ^ IS eas y *° show that be- 
low the incipient Stoner instability, J<J C , S± is domi- 
nated by the " white noise" term of Eq. Ijl9(l , while above 
this point it is the second (singular Stoner) term in Ijl9(l 
which dominates. 

Tunneling density of states. The conductance gx 
is related to the tunnelling DoS v through gx = 

— J dev(e)T(e) ( — <^T~) w ^ ere ^ F 1S ^ e Fermi distri- 
bution function at the contact and T is the golden rule 
dot-lead broadening. To obtain the TDoS from the GF, 
Eq. Q2U[l . we deform the contour of integration in accor- 
dance with 0. As a result, the TDoS is given by 



v{e) = 



+ itj)k,me iet dt. 

(23) 

where {...)k.m denotes a summation over all winding num- 
bers for Coulomb and longitudinal zero-modes [llj. We 
have computed the temperature and energy dependence 
of the TDoS for various values of e. These are depicted in 
Fig. 2. The energy dependent TDoS shows an intriguing 
non-monotonic behavior at energies comparable to the 
charging energy E c . This behavior, absent for J = 0, is 
due to the contribution of the second term in Eci. (|19|l . 
It is amplified in the vicinity of the Stoner point, and 
signals the effect of collective spin excitations (incipient 
ordered phase). 




E c /A = 10 



J I A =0, e = 



J/ A = 0.9, e = 0.8 







Energy 



FIG. 2: The spin-normalized tunneling density of states 
shown as function of energy. Insert: TDoS as function of 
temperature. 



Spin susceptibilities. The spin susceptibilities are de- 
fined through 



d 4 Z 



dK f dk» t dA.».dK i 



(24) 



The longitudinal susceptibility (x zz ) is not affected by 
the gauge bosons. By contrast, the transverse x H ac ~ 
quires the gauge factor (e 2ri ^) m ^z , where the average 
is performed with respect to the Gaussian fluctuations of 
$ z and, in principle, the winding numbers (cf. Ea.l|12|l). 
In practice, since T>J, only the m = winding should 
be taken into account; T>A allows us to evaluate the 
path integral in the Gaussian approximation. One finds 
to leading order in e 



X ZZ (r) 



Xo 



1 - Jxo 



X + -(r) = 



£Xoe 



Jr 



1 - tJxo 



(25) 



where xo — The above susceptibilities are given 

as function of t. To obtain the dynamic susceptibili- 
ties one needs to Fourier transform and then continue 
to real frequencies. \ zz (|2"5|) (lif diverges at the ther- 
modynamic Stoner Instability point, while x -1 remains 
finite. Notwithstanding, the static transverse susceptibil- 
ity is enhanced by the gauge fluctuations. The dynamic 
behavior (including relaxation processes) and the e cor- 
rections to will be discussed elsewhere. |l7j . 

Summarizing, we investigate influence of spin and 
charge zero-mode interactions on the TDoS and the sus- 
ceptibilities. Longitudinal spin fluctuations suppress the 
CB and the static longitudinal susceptibility is greatly 
enhanced near the Stoner instability. Transverse fluctua- 
tions generally tend to suppress the CB, but also contain 
a term which dominates the dynamics near the Stoner in- 
stability and enhances the CB. The transverse suscepti- 
bility will be enhanced as well. On a more technical level, 
Coulomb interaction is described in terms of Abelian 
(U(l)) gauge theory and lead to Gaussian gauge fac- 
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tor, the spin interaction, being a subject of non-Abelian 
(SU(2)) gauge gives rise to non-Gaussian gauge factors. 
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